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Abstract 

We have searched for lepton-flavor-violating r decays with a pseudoscalar meson (77, rj and 7T°) 
using a data sample of 401 fb" 1 collected with the Belle detector at the KEKB asymmetric-energy 
e + e~ collider. No evidence for these decays is found and we set the following upper limits on the 
branching fractions: B(t~ — > e~rj) < 9.2 x 10~ 8 , B(t~ — > /U _ r/) < 6.5 x 10~ 8 , B(t~ — > e~r/') < 
1.6xl0" 7 , B(t~ -> prrf) < 1.3xl0~ 7 , B(t~ -» e~vr ) < 8.0xl0~ 8 andi3(r~ fj,~ir°) < 1.2xl0" 7 
at the 90% confidence level. These results improve our previously published upper limits by factors 
from 2.3 to 6.3. 

PACS numbers: 11.30.Fs; 13.35.Dx; 14.60.Fg 
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INTRODUCTION 



Lepton flavor violation (LFV) appears in various extensions of the Standard Model (SM), 
e.g., supersymmetry (SUSY), leptoquark and many other models. In particular, r lepton- 
flavor- violating decays with a pseudoscalar meson (M° = 7),7)' and tt°) are discussed in 
models with Higgs-mediated LFV processes [l|, 0, y], heavy singlet Dirac neutrinos 
dimension-six effective fermionic operators that induce r — \x mixing j^] , R— parity violation 
in SUSY 0,0,1, type III two-Higgs-doublet models [s| and flavor changing Z' bosons jlj. 
Some of these models predict branching fractions which, for certain combinations of model 
parameters, can be as high as 10~ 6 ; this rate is already accessible at high-statistics B 
factory experiments. Previously, we obtained 90% confidence level (C.L.) upper limits for 
various r~ — > £~M° (where i~ = e~ or branching fractions using 154 fb" 1 of data; 
the results were in the range (1.5—10) x 10 -7 0. The BaBar collaboration has rece ntly 



used 339 fb _1 of data to obtain 90% C.L. upper limits in the range (1.1—2.4) x 10~ 7 [10 
Here we update our previous results using 401 fb _1 of data. These datasets are collected 
at the T(45") resonance and 60 MeV below it with the Belle detector at the KEKB e + e~ 



asymmetric-energy collider [11] 



The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL), all located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside the coil is instrumented to detect K® mesons and to identify muons 
(KLM). The detector is described in detail elsewhere 12| . 



Particle identification is very important for this measurement. We use particle identi- 
fication likelihood variables based on the ratio of the energy deposited in the ECL to the 
momentum measured in the SVD and CDC, the shower shape in the ECL, the particle range 
in the KLM, the hit information from the ACC, the dE/dx information in the CDC, and 
the particle time-of-flight from the TOF. For lepton identification, we form likelihood ratios 
V(e) [l3| and V(/a) [lH based on the electron and muon probabilities, respectively, which 
are determined by the responses of the appropriate subdetectors. 

In order to determine the event selection requirements, we use Monte Carlo (MC) sam- 
ples. The following MC programs have been used to generate background events: KO- 
RALB/TAUOLA 
events, KKMC E 



3 for t+t - , QQ (3 for BB and continuum, BHLUMI for Bhabha 
for e + e~ — > n + fi~ and AAFH [l9[ for two-photon processes. Signal MC 
is generated by KORALB/TAUOLA. Signal r decays are two-body and assumed to have a 
uniform angular distribution in the r lepton's rest frame. All kinematic variables are calcu- 
lated in the laboratory frame unless otherwise specified. In particular, variables calculated 
in the e + e~ center-of-mass (CM) frame are indicated by the superscript "CM". 



EVENT SELECTION 

We search for t + t~ events in which one r decays into a lepton and a pseudoscalar meson 
on the signal side, while the other r decays into one charged track with a sign opposite to 
that of the signal-side lepton and any number of additional photons and neutrinos on the 
tag side. Thus, the decay chain we reconstruct is: 



i 



{r~ -> £~ (e~ or yr) + M° (77,77' or 7r )} + {r+ -> (a track)+ + (n > 7) + X (missing) }[f]. 

We reconstruct a pseudoscalar meson in the following modes: 77 — > 77 and 7r + 7r~7r°(^ 77), 
77' — > p{— > 7r + 7r~)7 and r](— > 77)7r + 7r _ , and 7r° — > 77. While the tt° — > 77 and 77 — > 77 modes 
correspond to a 1-1 prong configuration, the other modes give a 3-1 prong configuration. 
All charged tracks and photons are required to be reconstructed within the fiducial volume 
defined by —0.866 < cos# < 0.956, where 9 is the polar angle with respect to the direction 
along the e + beam. We select charged tracks with momenta transverse to the e + beam 
direction, p t > 0.1 GeV/c while the photon energies must satisfy E 7 > 0.1 GeV (0.05 GeV) 
for the 1-1 prong (3-1 prong) configuration. 

Candidate r-pair events are thus required to have two and four tracks with a zero net 
charge for the 1-1 and 3-1 prong configurations, respectively. Event particles are separated 
into two hemispheres referred to as the signal and tag sides using the plane perpendicular to 
the thrust axis (2o| . The tag side contains a single track, the signal side contains one or three 
tracks. The track on the signal side is required to be identified as a lepton. The electron 
(muon) identification criteria are V(e) (V(p)) > 0.9 with p > 0.7 GeV/c. The efficiencies 
for electron and muon identification after these requirements are 92% and 88%, respectively. 
To reduce fake pseudoscalar meson candidates, we reject radiative photons from electrons 
on the signal side if cos 6^ > 0.99. 

The 7T° candidates are formed from pairs of photons that satisfy 0.115 GeV/c 2 < M 77 < 
0.152 GeV/c 2 , which corresponds to ±2.5 standard deviations (a) in terms of the mass 
resolution. We also require p w o > 0.1 GeV/c on the signal side. 

The 77 meson is reconstructed in the 77 (7r + 7r~7r°) decay modes. The mass window is 
chosen to be 0.515 GeV/c 2 (0.532 GeV/c 2 ) < m 77 (m^- w o) < 0.570 GeV/c 2 (0.562 GeV/c 2 ), 
which corresponds to —3. Oct and +2.5cr (±3cr ). To reduce background in 77 — > 77 channel, 
we reject those photons that form 7r° candidates in association with any other photon with 
E 1 > 0.05 GeV, within the vr° mass window, 0.10 GeV/c 2 < M 77 < 0.16 GeV/c 2 . 

The 77' meson is reconstructed in the fry and r]ir + iT~ decay modes. For the p — > 7t + tt~ 
selection, the mass window is chosen to be 0.550 GeV/c 2 < < 0.900 GeV/c 2 . We 
reconstruct 77' candidates using a p candidate and a photon on the signal side. The 77' 
mass window is chosen to be 0.930 GeV/c 2 < m pl < 0.970 GeV/c 2 , which corresponds to 
—3. Oct and +2. 5a. Furthermore, we veto photons from 7r° candidates in order to avoid fake 
77' candidates from 7r° — > 77. We remove events if a ir° candidate with invariant mass in 
the range 0.10 GeV/c 2 < M 77 < 0.16 GeV/c 2 is reconstructed using a photon from the 
77' candidate and another photon with E 1 > 0.05 GeV. Figure [TJ shows the p — > 7r + 7r~ 
and 77' — > fry mass distributions. The dominant backgrounds for this mode come from 
t~ — > h~ p°z/ r (+7r ) (where h~ = K~ or 7r~) with a photon from ir° decay, beam background 
or initial state radiation (ISR). As shown in Fig. [TJ there is no 77' peak either in data or in 
MC since decay modes with an 77' are very rare and are not included in the generic r decay 
model [ll[. We also reconstruct 77' candidates using an 77 candidate and two oppositely 
charged tracks consistent with being pions. We impose a 7 — > e + e~ conversion veto as 
-P(e) < 0.1 for both tracks in the 77' candidate. The 77 products from 77' — > r]7r + 7T~ decay 
are reconstructed using two photons. The 77' mass window is chosen to be 0.920 GeV/c 2 
< m^TT+Tr- < 0.980 GeV/c 2 , which corresponds to ±3. Oct. 



[|] Unless otherwise stated, charge conjugate decays are implied throughout this paper. 



5 



> 
o 



o 
o 



6000 



4000 



2000 



\^\t —> fj,r)'(—> py) 

• Data 
| | Backg round MC 




0.6 



0.8 1 
M„„ (GeV/c 2 ) 



> 
a) 
O 

to 
o 
o 



2500 
2000 
1500 
1000 
500 












< 


> 















0.95 1 
M„ (GeV/c) 



FIG. 1: The p — > tt + tt~ (left) and rf — > pj (right) mass distributions. While the signal MC 
(t~ — > p~T]') distribution is normalized arbitrarily, the data and background MC are normalized to 
the same luminosity. Selected regions are indicated by the arrows from the marked cut boundaries. 



To ensure that the missing particles are neutrinos rather than photons or charged particles 
that fall outside the detector acceptance, we impose additional requirements on the missing 
momentum vector p m i SS , which is calculated by subtracting the vector sum of the momenta 
of all tracks and photons from the sum of the e + and e~ beam momenta. We require that 
the magnitude of p m i SS be greater than 0.4 GeV/c, and that its direction point into the 
fiducial volume of the detector. Since neutrinos are normally emitted only on the tag side, 
the direction of p m i SS should lie within the tag side of the event. The cosine of the opening 
angle between p m i SS and the thrust axis (on the signal side) in the CM system, cos6^ s _ thrust , 
is therefore required to be less than —0.55. 

For the 1-1 prong configuration, to suppress fake r\ (ir°) candidates arising from beam 
background and ISR, we require that the higher and lower energy photons (E 7l and E l2 ) 
in an t] (ir°) candidate satisfy the requirement E y i > 0.6 (0.9) GeV and E l2 > 0.25 (0.2) 
GeV, respectively, as shown for the r~ — > fi~TT° mode in Figs. [2] (a) and (b). In order to 
suppress background from qq (q = u, d, s, c) continuum events, we require that the number 
of extra photon candidates on the signal and tag side (n^ IG and n^ AG ) be n^ IG < 1 and 
n TAG < 2 ; respectively. To reduce background from Bhabha and events, we require 

the momentum of the lepton and that of the tag-side charged particle in the CM system to 
be less than 4.5 GeV/c. Furthermore, we require the momentum of a lepton to be greater 
than 1.5 GeV/c for the r~ — > £~n° mode (listed in Fig. [2](c)). This condition is not imposed 
for the r~ — > £~r] mode, in which the average lepton momentum as well as the background 
level are lower. 

The reconstructed mass on the tag side using a track (with a pion mass hypothesis) and 
photons, m tag , is required to be less than 1.777 GeV/c 2 . The total visible energy in the CM 
frame, E°-£, is denned as the sum of the energies of the r\ candidate, the lepton, the tag-side 
track (with a pion mass hypothesis) and all photon candidates. We require E™ to satisfy 
the condition shown in Table [B To reduce background from fi + fi~, two-photon and Bhabha 
events, we add the veto condition E™ > 8.5 GeV for the muon (electron) mode if the track 
on the tag side is a muon (electron). The cosine of the opening angle between the lepton 
and the M° in the CM system, cos0™ M o, is required to lie in the range shown in Table [J 



6 



> 

CD 

O 



o 
\ 

c 

CD 
> 
V 

H— 

o 
-1J 



o 
> 

1) 
o 



c_ 

CD 
> 
CD 

M — 

o 
a; 



400 
300 
200 
100 




(a) O^" 

I" 1 J • Data 

| | Background MC 

> 




> 

'CD 

° 500 
o 

o 400 

\ 

r 300 
c 

> 200 
o 100 




(GeV) 



<D 

E 









CM 

92500 
6 

>2000 

=t=> 

1 1500 

CD 



4 6 

p u (GeV/c) 



0.4 0.6 



2 

E )2 (GeV) 



- (d) 






<: 1 











0.J 
cos 1^1. 



FIG. 2: Kinematic distributions used in the event selection: (a) higher energy and (b) lower 
energy of a photon from the 7r° candidate (E^i and E^)] (c) momentum of a muon (p^); (d) the 



cosine of the opening angle between the muon and 7r in the CM frame (cos# 



CM 



While the 



signal MC (t~ — » /i _ vr°) distribution is normalized arbitrarily, the data and background MC are 
normalized to the same luminosity. Selected regions are indicated by the arrows from the marked 
cut boundaries. 



(also shown in Fig.[2](d) for the r~ — > /i~7r° mode). For all kinematic distributions in Fig. [21 
reasonable agreement between the data and background MC is observed. 

The correlation between the momentum of the track on the tag side, p™, and the cosine of 
the opening angle between the thrust and missing particle, cos ^miss-thrust 111 ^ ne CM system is 
employed to further suppress backgrounds from generic t + t~ and fi + [J>~ events via the follow- 
ing requirements: pgf > 1.1 log(cos^M_ thrust + o.92) + 5.5, andp™ < 5cos^ s _ thrust + 7.8 

where is in GeV/c (see Fig. [3]). Finally, we require the following relation between the 
missing momentum p m i SS and missing mass squared m miss to further suppress background 
from generic r + r~ and continuum background. In signal events, two neutrinos are included 
if the r decay on the tag side is a leptonic decay, while one neutrino is included if the r 
decay on the tag side is a hadronic decay. Therefore, we separate events into two classes 
according to the track on the tag side: leptonic or hadronic, and apply the requirements 



7 



> 



CL 




> 

O 



2l 



4 

3 
/ 
1 





(b) Data 















ii£^* .*.*> 


•* .' ;} •• ■' ■ . ■' ' 


~! »* .».*«# 
:»•• - *••../>• 


' " ■ • 7 *• 


Pp.;.' 




eK'-'r*";'' 
PsT. ' "•• 


I" " ' 1 



0.8 -0.6 

COS 1^ m j S s-thrust 



-1 -0.8 -0.6 

COS ^miss-thrust 



(c) generic tt 




-1 -0.8 -0.6 

COS miss-thrust 



FIG. 3: Scatter-plots for (a) signal MC (r~ - 
and (d) MC events in the — cos # : 

curves. 



(d) /llja 




-1 -0.8 -0.6 



CiJS ^miss-thrust 

fj,~r)(— > 77)), (b) data, (c) generic t + t~ MC events 
is-thrust plane. Selected regions lie between the two 





771CM 


cos ^^[^0 


T — > £ T] 

t~ r?/ 


0.529 GeV < E% <10.0 GeV 
0.529 GeV < <11.0 GeV 
0.529 GeV < £^ < 10 - GeV 


0.50 < cos 9f^ v < 0.85 

0.50 < cos 0f™ v , 
0.50 < cos 6™ w0 < 0.80 



TABLE I: The selection criteria for the total visible energy in the CM frame (E™) and the cosine 
of the opening angle between the lepton and the M° in the CM system (cos^r^o). 



shown in Table [TT] (see also Fig. |4j. 

For the 3-1 prong configuration, we impose similar requirements m tag < 1.777 GeV/c 2 , 
n siG < £>cm anc j CO s9™ M0 (Table[I|. For the 7/ — > fry mode, we require the photon energy 
to be greater than 0.25 GeV for the barrel and 0.40 GeV for the forward region to suppress 
fake candidates from beam background and ISR. The cut on the n^ AG and jj^-cos ^miss-thrust 
correlation is not applied for the 3-1 configuration modes, in which the background level is 
lower than in the 1-1 configuration modes. We apply requirements on p m i SS and ?n miss similar 
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FIG. 4: Scatter-plots of p m i ss — remiss' f° r leptonic and hadronic tags: (a) shows the signal MC 
(r~ — > [i~rf) distribution with arbitrary normalization for a leptonic tag while (c) shows the data 
distribution for a leptonic tag; (b) shows the signal MC (r~ — > n~rj) distribution with arbitrary 
normalization for a hadronic tag while Fig. (d) shows the data distribution for a hadronic tag. 
Selected regions are indicated by lines. 



to those for the 1-1 prong configuration (see Table 1171) . 



SIGNAL REGION AND BACKGROUND ESTIMATION 

Signal candidates are examined in the two-dimensional plots of the £~M° invariant mass, 
M inv , and the difference of their energy from the beam energy in the CM system, AE. A 
signal event should have M- mv close to the r-lepton mass and AE close to zero. For all modes, 
the Mj nv and AE resolutions are parameterized from the MC distributions with asymmetric 
Gaussian shapes to account for initial state radiation and ECL energy leakage for photons. 
The resolutions in M- mv and AE are given in Table IIHI 

To evaluate the branching fractions, we use elliptical signal regions, which contain 90% 
of the MC signal events satisfying all cuts. These signal regions are shown in Fig. and 
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TABLE II: The selection criteria for the missing momentum (p m ; ss ) and missing mass squared 
( m miss) correlations, p m i ss is in GeV/c and m miss is in (GeV/c 2 ) 2 . 
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23.1 


35.6 
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TABLE III: Summary of M inv (MeV/c 2 ) and AE resolutions (MeV) 



the corresponding signal efficiencies are given in Table IIVI We blind the signal region so as 
not to bias our choice of selection criteria. Figures [5] and [6] show scatter-plots for data and 
signal MC samples distributed over ±10<r in the M mv — AE plane. Most of the surviving 
background events in r — > £ir° modes come from r _ — > 7r~7r°z/ r , where the tt~ is misidentified 
as a lepton. The remaining backgrounds in the r~ — > n~r](— > 77) mode are from r decay 
including a real 77 meson or combinations of a fake lepton and a fake 77 meson formed by 
7's from 7T° decay, ISR or beam background. As there are few remaining MC background 
events in the signal ellipse, we estimate the background contribution using the M inv sideband 
regions. Extrapolation to the signal region assumes that the background distribution is flat 
along the M- mv axis. We then estimate the expected number of the background events in the 
signal region for each mode using the number of data events observed in the sideband region 
inside the horizontal lines but excluding the signal region as shown in Fig. \5\ and O The 
numbers of background events in the 90% elliptical signal region are also shown in Table HVl 
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TABLE IV: Results of the final event selection for the individual modes: Bj^o is the branching 
fraction for the M° decay; bo and s are the number of expected background and observed events 
in the signal region, respectively; "Total sys." means the total systematic uncertainty; S90 is the 
upper limit on the number of signal events including systematic uncertainties. 



Systematic uncertainties for M° reconstruction are estimated to be 3.0%, 4.0%, 4.0%, 
5.0% and 3.0% for rj —>■ 77, 77 — > 7r + 7r _ 7r°, rj' — > pi, rj — > 7]tt + 7t~ and 7r° — > 77, respectively. 
Furthermore, the uncertainties due to the branching fractions of the M° meson are 0.7%, 
1.8%, 3.1% and 3.1% for 77 — * 77, 77 — > 7r + 7r~7r°, 77' — > p7 and 77' — > 777r + 7r~, respectively 
211 ] . For the n° veto we assign a 5.5% uncertainty for the 77 — > 77 mode while a 2.8% 



uncertainty is assigned to the 77' — > /ry mode. The uncertainties in the trigger (0.5—1.0%), 
tracking for lepton on the signal side and track on the tag side (1.0% per each track), lepton 
identification (2.0%), MC statistics (1.0—1.5%) and luminosity (1.4%) are also considered. 
All these uncertainties are added in quadrature, and the total systematic uncertainties are 
shown in Table IIV1 

While the angular distribution of signal r decays is initially assumed to be uniform in 
this analysis, it is sensitive to the lepton-flavor- violating interaction structure 22| . The spin 



correlation between the r lepton on the signal and that on the tag side must be considered. 
A possible nonuniformity is taken into account by comparing the uniform case with MC's 
assuming V — A and V + A interactions, which result in the maximum possible variations. 
No statistically significant difference in the M inv - AE distribution or the efficiencies is found 
compared to the case of the uniform distribution. Therefore, systematic uncertainties due 
to these effects are neglected in the upper limit evaluation. 

We open the blind and find only one event in r — > fnr°(^ 77). In other modes, no 
events are found in the blinded region. Since no statistically significant excess of data over 
the expected background in the signal region is observed, we set upper limits for branching 
fractions. The upper limit on the number of signal events at the 90% CL. sg including 



systematic uncertainty is obtained with the use of the Feldman- Cousins method [23j calcu- 



lated by the POLE program without conditioning [24J. The upper limit on the branching 
fraction (£>) is then given by 
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TABLE V: Summary of upper limits on B at 90% C.L. 



where Bm° is taken from Ref. 
in 401 fb" 1 of data. We obtain N, 



2l| and N TT = 357.7 x 10 6 is the number of r + r pairs produced 



using cr T 



0.892 ± 0.002 nb, the e+e" 



t + t cross 



The combined upper limits for 



section at the T(4S") resonance calculated by KKMC [18 
the r) and rj' modes are obtained by summing e£>, the observed and expected background 
events of each subdecay, and systematic uncertainties are estimated by first summing all 
correlated terms linearly and then adding quadratically the uncorrelated terms. The upper 
limits for the branching fractions B(t~ — > £~M°) are in the range (6.5 — 16) x 10~ 8 at the 
90% confidence level. A summary of the upper limits is given in Table |V] These results 
improve our previously published upper limits [o| by factors of 2.3—6.3. They also improve 
upon the recent BaBar results [Toj by factors of ~ 1.5. 



DISCUSSION 



The branching fraction for the r~ — > fi~~rj mode may be enhanced by Higgs-mediated 
LFV if large mixing between left-hand scalar muons and scalar taus in the corresponding 
SUSY model occurs [lj. This can be written as 



B{t~ -> iTrj) = 8.4 x 10' 



tan/3\ 6 ZlOOGeV/c 2 



60 



m A 



(2) 



where rrtA is the pseudoscalar Higgs mass and tan j3 is the ratio of the vacuum expectation 
values of the neutral Higgs fields coupled to up-type and down-type fermions. From our 
upper limit on the branching fraction for the decay r _ — > fi~r], some region of and 
tan j3 parameters can be excluded. Figure [7] shows the excluded region in the — tan j3 



plane. It also shows the constraints at a 95% C.L. from the CDF [25J, D0 [26J and LEP2 
experiments 27J. The excluded regions from these experiments are shown with the Higgs 
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FIG. 5: Scatter-plots of data in the M- mv - AE plane: 
±10(7 area for the r~ — ► /z~r/(— ► 77), r~ — » e~rj{-^ 77), r" 



1.75 



^(Gev/cf 



a), (b), (c) and (d) correspond to the 



A* »7(- 



7T 7T 7T 



and r 



7T 7T 7T 



0> 



modes, respectively. The filled boxes show the MC signal distribution with arbitrary 
normalization. The elliptical signal region shown by the solid curve is used for evaluating the 
signal yield. The region between the horizontal lines excluding the signal region is used to estimate 
the expected background in the elliptical region. 



mass parameter \i > in the maximum stop- mixing scenario 28]. We note that their 
theoretical assumptions are somewhat different with from ours, and thus these regions are 
for illustrative purposes only. 

The improved sensitivity to rare r lepton decay achieved in this work can also be used to 
constrain the parameters of other models, e.g., those with the heavy Dirac neutrinos jij. In 
this model, the expected branching fractions of various LFV decays are evaluated in terms 
of combinations of the model parameters. These parameters, denoted y Te and y Tfl for r 
decay involving an electron and a muon, respectively, can vary from to 1. We obtain the 
following upper limits: y re < 0.17 and y TfJi < 0.47 at 90% CL. from our r~ — > £~tt° results. 



SUMMARY 

We have searched for lepton-flavor-violating r decays with a pseudoscalar meson (77, 
rj' and it ) using 401 fb _1 of data. No signal is found and we set the following upper 
limits on branching fractions: B(r~ — > e~t]) < 9.2 x 10~ 8 , B(t~ — > ^ff) < 6.5 x 10 -8 , 
B(t~ -> e-rf) < 1.6 x 10- 7 , B(r~ -> fTrf) < 1.3 x 10" 7 B(r~ -> e'n°) < 8.0 x 10~ 8 
and B(t~ /i~7r°) < 1.2 x 10 -7 at the 90% confidence level, respectively. These results 
improve upon our previously published upper limits by factors from 2.3 to 6.3. They are 
also somewhat better than the recent results from BaBar 1(| with the single exception of 
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FIG. 6: Scatter-plots of data in the M- im - AE plane: (e), (f), (g), (h), (i) and (j) correspond to the 
±10<7 area for the r~ — > /i - r/(— > pj), r~ — > e~r/(— > pj), t~ — > n~r]'(—> 7]ir + ir~), r~ — > e _ r/(— > 
t]7t + tt~), r~ — > /i~7r°(^ 77) and r~ — > e~7r°(— > 77) modes, respectively. The data are indicated 
by the solid circles. The filled boxes show the MC signal distribution with arbitrary normalization. 
The elliptical signal region shown by the solid curve is used for evaluating the signal yield. The 
region between the horizontal lines excluding the signal region is used to estimate the expected 
background in the elliptical region. 



the limit for the r — > /x n mode, and are the most stringent limits on these modes to 
date. These limits help to constrain new physics scenarios beyond the Standard Model. 
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